{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "injured-accident",
   "metadata": {},
   "source": [
    "# PySCHISM example_1.py\n",
    "### Example file for basic model deployment.\n",
    "---"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "lucky-doctor",
   "metadata": {},
   "outputs": [],
   "source": [
    "from pyschism.mesh import Hgrid"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "greatest-curve",
   "metadata": {},
   "outputs": [],
   "source": [
    "hgrid = Hgrid.open('../../../NOAA_project/hgrid_atlantic.grd', crs='epsg:4326')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "faced-ecuador",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAOEAAAD4CAYAAAAaaoEpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABUbElEQVR4nO29f5hcV30f/Dkzs7Oj0Wg0Hq2WlbysZbEWQghFGCEcA45rHNclBBzq8gCvQ3hLAvQJb0IpaSDJ2/I0TUpSUpI+NORxoAGatoY6NHF5Er8EiOu6fg2vrChCCKMIRSyyWK+XZb1eRqvR7Jz3j+/97Pnes+fOj90rjWTdz/PMMzN37o9z75zP+f4832OstciQIcPgkBt0AzJkuNqRkTBDhgEjI2GGDANGRsIMGQaMjIQZMgwYhUFdeGRkxO7YsWNQl89wFeHxx+X9ZS9L+5yPz1prt673XGsioTEmD+AQgCetta8zxnwKwI8BeCba5e3W2iOdzrFjxw4cOnRoLZfPkKErjIl/X08kTp+rUADabSCXA9pt8521n1Wdc43H/SKAbwKoqm2/ZK29f/1NypBh7fDJ12l7L8TUx+Vy8ipErFla6r99IfRtExpjxgH8BIBPpNOEDBn6hzHhF+DIkgRruxNQn0+j3QZaLXmlhbU4Zn4XwD8H0Pa2/4Yx5qgx5qPGmOHQgcaYdxpjDhljDj399NNruHSGDKtB6VQoUE2UF+AIp19J8Mmsz10qAcVifHtq7e9nZ2PM6wDMWGsf9376IIDdAF4OoA7gl0PHW2vvtdYesNYe2Lp13fZshqsQPklIPkCkk0++0PH6PJ2IVyzGic3z0yZMC/3ahK8E8HpjzGsBlABUjTF/bK29J/r9vDHmjwC8P70mZrjaoAlhbbKdBzjSUT1MknT6HKHzUYX1VVmeX0tXfb000BefrbUftNaOW2t3AHgzgK9Ya+8xxmwDAGOMAXAXgGPpNTHD1QSfICEJxVexKNu6EdCHdrAUi07V5Pn8fUk+/dsgJWES/rMxZisAA+AIgHendN4MVxGSJJ4vpdq+N6KP81K91N8JLfVC7/ysiZ8G1kxCa+1DAB6KPt+WUnsyXEXopGaSKFo99N99FdGY1dLQ95hqEtLLWSyGiZdE9n4HgW4YWMZMhgw+fHXPt9EYKNcE0YQEHBFDamxImlKqdSJdUhvTImNGwgwDgS+hQjZWiIDcDiSTICleqD2bJKZPav/6IXuQ35vN7vfZCzISZrjk0CTRNpkmGn/3Cder9KHzRttunY5NsjV5/UpFXoUCUK8D1Srw0EO9taVrW9M5TYYMvSGJgASJwN98Uvrw7US9vV91UZ+Lx9brwMiIkK5add/L5YyEGZ4DoGpICehLrU6ka7edpNMqY4iM/ucQcYtFIZa2HUm68XFg+3YhX70O1GqXZ4giQ4Y1QztLfEJ1ssm0feerq6H9q1VHtHJ5tc1ZKgnBGDusVIDRUUfAatUNGtUqUkNGwgyXDH5IImT/JYUItNTk/txHS1JfgtI23L4dGBuLk69UkhcJWCgI8UjWctlJv0pFvrNNGQkzXPHQZGu14jZgSBVtt8Ub6dt5SY6dXE6k2o4dok6Wy0IkZsUwW4ZSr1SS7ZWKkK5cdpKxWl2dUbOpvJzas8hImOGSgYF0XyI2m+FwAsMAnRwsvrSkqjk5KQSs10WlpAQkmYD4d5KxWhXikYSUiLkckG9fcKPGQkrxCWQkzHAJQNIlJWNTyvkpZCECJjlEtMNmxw7g5pvjEo22oJ4dob9zH4YiymXZnm+eAxaX3MlJwrSChMhImOESolOaGu25UGiiG4pFsd0AJwXHxhwJSSitgmoy8p37lcvAUPu8kK/ZjM/i5ecUc9cyEmZIBZ0I1g+6eTk1+Bsl2I4dQr5du5xDpVqNO2C0bed/L5UilbPZlNoVS0vxnDYtCVPM4M5ImGHN6Id43VLN9G+9xOB06hmnIo2OAnv3CgHHxlyAnSqplnw8x1DBOnI1IsnHl09AIJOEGS4P9Cv1kubhrbUfawJWKsDEhEi/ffuA3bvj5NMkzOdsnFitNrDk2Xl86ewBP26SkTDDINErAXuVZr2UigjNrigWhXgTEy4MsXOnxANrNaeCrth3S62waqlJpQmYNLcpZQIC6dUdvR7AfQC2AHgcwE9ba9NzH2W4ItApzcxHp8wW/U5HjZ4HWKvJa3JS1M9aTVTR8XEXjiiVgOHCMtBohG07NkITi6omJWAoe6DTtIs1Iq26o78F4KPW2vuMMX8A4B0APp5C+zJcZgiFGbqFDfxt3aSfnv9HhwrTznbsEJuvUBDJNzEhJGRWS0z6LTTCDhaik1TUjQwRd5AkVHVHfwPA+6K6MrcBeGu0y6cBfAgZCZ/z6Jd8STMedH9nzK5QELLpbJd2W4i4fbtIPCZT0wNKEg7hgiNfJxLy4nwPSUu9X8rkI9YiCX8XUnd0U/R9C4B5ay19tmcAXBs60BjzTgDvBICJiYk1XDrD5YReJJp2oDDzBBCy0QHp78eslZ07RfJRJWVmzciIkLBadcQjCTeULNCIQgy+oyUkCfmepG767/7nFNAXCXXdUWPMrf1ezFp7L4B7AeDAgQPZOt2XGXRmS7d9QvAlXbEooYKdO13+ZbXqAuMLC0Ke+XknhCoVIRlnLTCAzhADpaJ2vlBy5nM2bgOGCOg7Vnp1wISKm6aEddcdBfB7AGrGmEIkDccBPJlqKzNcdCSVGtSE7NUW5FQfEnB0VD7rRGnA9fWJidVxP+3d1CUJGXZgelm+eQ5otiTUADjiURImSbkktdTf5je60Ug1ZQ3ok4TW2g9Cqm0jkoTvt9b+H8aY/wbgboiH9GcA/FmqrcxwyZBUrcxHkleTTpTRUQkh8J0SsNP8Pz+LxSegTrLeVLHA4qLTiZntTfJpO5AX6kY8v1Ga7YCcU8cPU0JaccJfBnCfMeZfA/hrAJ9M6bwZBoROla+TCFgqiVSbnHRq5cREPCG6kzOHaiVfWiJqIm4qLwMLi/E5UIAjWkgCJt1AN48p9/cnKV4OaWte3dFTAA6m06QMlxprST/zwdzNAwecDcgyEP6MBcD1b0pHPbFWk1DbhIVC5PlcWHRSSS8WoSWgjv3xgppsOvG0UBAD1Q9NLC2tHjmSApzrQJYxkyER3daBAFz8bvdu4JZb3OwFloLQ2iKJpAmopxX5sxq0JMy3L8RVTZ1qBsS/axKGPJ25nIwU9TqWC8Pi0JmaAubmZJ/Fxfh56EmiQTo9LfukhIyEVzFC8/xCVawJf/AvFoVw+/e7+Xv1uiMbj9HBdy3pfPIxHU3P8RsuLIfDDVrSdQsv8LPeXijgydlhzM4CtZrBdQyZLS7GyVcqAWNjeLa4Ba0WcE39gmyfnl738ycyEmYIeka72YSFgjhdDh4Ux0utJjYgJRf347uWhnqCrX75jpl8K5rTp8MNIVUTSHbChDJiVA4cSb/cNshrzxFHkcgg5eWXq0PI6+n5KSAj4VWMTuGHbmroyIjkbW7fHg/C6wRr/zuD8dr28/dbsfuWvClF/vSiUFXfJCno7wcA09O4dqSFa3fWRNKenI6fQ4vyRgNbKjmgUgBm5kRtHR3t/ID6QEbCDD1Bk2xszE0bqtXiDkqfhEC8fATDDiRgPqdGAjpXfOIl2Xm9ZLWEPJ2tFnD2rKiUzIfTOrM+z1xEOj3tfvt2LFevSeW5AhkJM/QBOmE4hYgeUFak9tf7Y99mqI1EHC4sR9IsQKIk6dfJ9tNIsgOTJjU2Gm6bJqH/ajZdbly5jPn5tJ5qRsIMPUCnoe3bJ6EInTamY3haAtKpyJh3sRipms2ERGlNwpD62Y2APEeo8Tr0kLSv35ZQilv0Kpe7PbXekZEwQ09xQiZOMw2NEpDk8gso+RXMikVgKLcsBPSJpeGHGTqpof2gU6wv6Zw6vUe/5uexoZxVW8swAJCAJKG28zQJNfn4mxAwkFDt23b+lCPfs8ltPrrNKO5EXD/z3M+Q0YPGzIyosCmKwoyEGTp6SQHnLKxU4o4XTT6/olks0I7lzqolL+K/+xLLX1AwaYJiKDk1CZ2m9uvPeiCYn8+C9RkuLeiQGR11NiCJqBdP0eUDV4orYbm7aplgd61cPNQg/Xu3fTrdWOi7T8AQqdeiEicgI2GGGJKC9HohFQbdQy+d9zmUW47beCHyJUlAPxao7bckhMjYyYvqE5nStRPB9X4pISNhhp5BKafTz3TsT9eDyeWwWvXslmIW2q9X+CGIJM9nt3NqIvquXr1PJgkzXExoaci+PTYWr1ydNNWIv+XbFzpXONPb/FQ0vvvThfrp/D4RQ1JX36B/nJ/82mn60zrRb3mLEoCHAQxHx95vrf2XxphPAfgxAM9Eu77dWnsktVZmGDhaLRfzY6K2TrqOOWF6Da5rsoWI2MlpA3T3iIY++8cm2YE6odUnYS/2Zh/oVxKeB3CbtXbRGDME4BFjzF9Ev/2Stfb+1FqW4bIAtbOdO11peSZo6/zPlYRrv0x8SAr6xEty3PAY/z2JUJ2QRE5+96d7+J+7nWMd6Le8hQVA3+xQ9MoKNj1HkBS037fPzRVkiUEdA1xJuvZXMArF/zqRT88NTFIfNbQ6q2c+EElE6aaOJknAi4S+r2CMyRtjjgCYAfCX1tqvRj/9hjHmqDHmo8aY4YRj32mMOWSMOfT000+vvdUZUkeIgAxN/OzPuiK727c7AtIu3FhWVc503qfvbPFT0vQ8QX/OYLfJuZ0I3c2D6hPMV0m1OqrV0ouEvs9srV221u6HVFU7aIzZCyn+tBvAywHUITVnQsfea609YK09sHXr1rW3OsO6YUz85YOzJX77t12yNos1sU+y9CAWFoSESbFALQVZjLfRiL86rYbU7ZVkVyYV8w0Rkdv17/5Ex4tExPXUmJk3xvwVgDuttR+JNp83xvwRgPen0roMqaOXchWABObf/W4hop4hwcwZmmYrU5FyubhKqZ0Xer2/JEnnOzt6DcBrTyYn7FI9Dd2Yv13blZ2IehHqjRL9eke3ArgQEXADgB8H8FvGmG3W2u9FJfHvAnAs/aZmWC96TdQuFoG3vhXYs0ckHbNh+PmaapSGVizKdCQtUdrteF1OfvdJmFSEV3/vZpNpArLx/s0kkUtf06+x4QfrfYmasq3YryTcBuDT0apMOQCfs9Z+wRjzlYigBsARAO9OrYUZLir8vhSVVFlJwgakbzIssQHngCbiBZe0zec7WaiChurDJEmWTp28kx0X8nJyn5Anlftu3+5KhLfbMol3aUn28b2xXnmMNNCvd/QogJcGtt+WWosypIZeVU9+Zpn6yUmpHcPQAyflbmh83x3kS4eQmhn6nBQ37NSpQ4RLOs4noq9u+lKwXgdGRvD9xgY0o8kRm0cLQsRGw+2vg/f6+BSQZcw8R9FvLVFWTqtWhYCsmlaruVAEylXpmKzRSfiz4H0yhjyXvXRi34PZzZniezO1OhoiOdXZQmFFsBeLiEs6TWJdg9TP5lkHMhI+B9Gr7Qe4MAPDDrt3y4t9TptKaLWcN1NDq5whyddL6EA3LES0TiT0PZhaEhI+sdiWxUVgcRH1+jXOxFxYkOlKSTHEFAkIZCR8zqEfCagXXCmXRepNTLg0NE7cLRaBYZx3nZNrlGnvZ0gV7ZSyFkIv5PPtMT+UoInGkSR0bm3Lzs9jeHHRnYPFf0OLZwDhvNZ1ICPhcwi9LlvmJ2BzFny97koYUjpypgQaS24yq+7QOvbXaykKXz3s1c4LkQ1YLQkp0nhTpZK0T6vR2tPZbIp01+31g/ShRIGUkJHwOYJui7eEiuySiHQM7tzpZs9T61oxg+h8YeejqtctuwXoHIoINdb/3GnfJBW1XAYmJvDXRwxmZ4FKZQh7927CpoUn3T6hdur780MW/Ur2HpGR8ApGL97PkL+Cn+n11MtOa9Kyn+Vb50WS8MBQnI9IyvvspdMmOVD8cyd5OrmtXAZ27MCnP2PwiU/IQLN3rzT95du7tCFJ7dXXTxkZCa9A9Bp6YF/SzhWaSVQ5q1WxA3fudLPnqb2125FJlfMcK1oV8+Nva5UUPgE7eTT9c2tpFS0P/Ou/mcfHPiYaaLns1rhfWac7Cb5UBcJzHTNJmCEJWpKRbNzulyOs1SQwPzERT8iemZE+VirJe3XnRgy1ZzsTMST99O/9oJNE1NcJ7TM5if99eAMefthJ+9tuA977XuAF288B84EZHoR+eH48NDTfMSVkJLzC0IvzhVpjteq2+zPf+eIiLiQh/RTz80LSlb7m13vRNlUaHTNJ6mmxzvdOUrPVwrFjQr79+4G77pK458iIOtYvuRGCb/+FEsRTQkbC5wB8RyE9m37oTHNEhyV0qUJ69LXdONQ6t/pCvsOCSIOISWGK0Hd//5kZ3HLLphXpvmcPsLVyTtq10IgTSMcStYT1pbwfdskk4dWLXmxBqpm6P/lefRJsbCxexlBP0uX3chluVoJmc6+qWb+dtVOcsNP+fG808KLit/CiPVyRtAA0FNFCQXuSTktIfa96BeAsRJHBhy8FK5X49pCHVM8H9Av1+ktWLy0BqFfcnEHC75AXyXO4gtC5Q6op99XJAkwu0Pmf2suqz6+nVhHdQjDrREbCKwTdpCAJSGeKDkcAzi4kKfVKSVRfGZwHnCTM5dQCmmRoqSQdc2FBgvcXQToEs1/4mb+HCOjbrkkOo1CbfXXUv6ZGRsIMQLz/6XUg+BuJp/uUJqG/ngSD9tyXmTT55jk5CVnabouILJdltsHs7MWRgp08rv4DCNml/vH+PEedfuY7mnK55PS0lO81I2Ef0Gu8h74PClqy+V52SkUfnAyhyxeSZwza53LApooFFltOrGpp2G7LiUql1bMlujW4X3SKGybtnxTP1GqqnoWv1U3thOHxIXKngL5I2KHu6PUA7gOwBcDjAH7aWpve2lGXGL0kQet9/P0vNSn9oryE3+/4mb+dPevIWi6vTK3Dlsp5oFDAMvJxJ4WvHtKdWi6HCzytF91ihZpAuo1sA6W2LzFDcc1Q6EKn6ekHOmCbMKnu6PsAfNRae58x5g8AvAPAx1Nr5SVEvwTs9vulIKQOL4Ti2LqvUpNcXBRNkgRlOcOVEzabyFMl07PktarWbjvPTi7n0mzY+TupkutFS0ln/4Z1W301wLcTQ55S30O6lrmQfaAvEnaoO3obgLdG2z8N4EO4AknYzzQgjSTThOfstvRYL0T1F2rRAiA25w+r+5A+hr8tRZMiGNTXWtgPFvK4plZyU3rootedXnsMy2XX4RkTWY+zppO6GiKbvq6Gbwf6x/sBVP1d6/U8Vz/zIvtA3zZhVF/mcQCTAP4DgG8DmLfW8mmfAXBtwrHvBPBOAJiYmFhLewPnXL0tqdOvRyqFHHXcrjt3yKG2VnJ3gt9PQnZfp/YC0j8XF4HxcQlqs7ZMqQRXllBPVdIXosqnl2vyO6wfBugHnWKDoZsEVkvHJK9m6NxanU2auOtLy5TQNwmttcsA9htjagD+O6TeaK/H3gvgXgA4cODAminRjzoY2h4iY69E8eNy9Ems57/xpWWv8LWjJOnnv+uBoliUNehfffC8m8i6uCheGyBuG+mJroDLAg/NvdNFoNaDUAywm4NG69/6d71PaF6i/s7Bw4+BXgTVes3eUVV39EcB1IwxhUgajgN4Mq0G+khbqvRTCkKDahxDZizOdTGJmLR2oK9xhXwp+p2CQkcb0GwCJ0/GpYCv82pJwCAj02y0N5EE5LG9kjFpnyQ7LkQ0rRKE1BI/DJEkbf2SHaHrp4R+vaPBuqMA/grA3RAP6c8A+LNUW4l0yddvESRfUykWxYuoZ8XwPw0VFeu3Xf1KRfaV0OCu74Pv7KesK7O0BGfX6RugqKf7Vbvu/Y7sT9HwxTM/hzpyt87tP0St5vo2nZ9iFyKqHzj14Zfm9501KaNfSZhUd/Q4gPuMMf8awF8D+GSajbwYNlWvYL8aGXHVyOp1+X9mZ0VrC9lknXwS3YjZSSqGpCGlGXnk9y1eT2uMDOK3WtFsi2ZTjMKZGdlBOzpCDgvf8aJtQD4QTQi/wfohhH7X8A1yfvdDFL4K4BPRTz0LPSQtzXWtVH/QSZGM/XpHjyJcd/QUgINpNYq4lOQL2enlspCvUpFEZ04MLRTEvb+46NIpQ3Z9SBvivmmZFtqr2Wq5aIHfpiRi1mrRNJ+QZ0cbu36jaSOyQBLgFoUJpeh0uuFuaia3+VNCfLutExH9++F19Tm43S/feDmR8FJhLU6SXtCpH1Cb4v9QKgnpONOgVnOSkM5CPeMg5BzR5lWaRPSlIfsQVVKmnvmCKQQG6VcKIflpXNS3/dhfqyUE1PERf6ZBL8TjtTqpovxd1/3UN+6rmPp3LfFC+/oSUL/8so1XAwl7mbC6HiSdQ//HzKdkUnO16gjItRg4W10TlnE3fR0dwgJ6J916QilaqCTZiJofExPR7zoUEWooXamhuJl/Yv/Va8N7+ZN5Y9o+TTof9w+1RxPUJ7P/0r/rNqSEy4aE3aqF9ftbP/tTUjADi1Xy+JnEo4+Ci2SyRCBzmKennTCZnY1rZiEipqmWaimo+7MmYej9kUfkfnbccw2G26dch9NTLviQNAm1iup3bA1/W+i7fhhrGW2THDdsZ7djQ8TrRsLnmiT0CRiyz3x0k2pJv/nv7GskYLXqZg/oqT76RSKOjTnhMTMjJFxclJzMYlEIOTcn1wr1u27Po5NETApXsCyodlIylBfSpJaWpI2tFjBcq7mgvM9erWNrPTdJgvRKPv09JLJ9lVETIklih96TrqnPmZSErvcJVWFbJwZOwhABk0iotyeRTfcVfYzuN9xHl3XwZ5NT8vmTXEP7N5vi3Ni1S/rw9DRw5AjwxBNu2l3If8D2rVca+v2Ic1i5mIvffl8aV6tyzMZqVUS4Tz4djA81PKSC9kIU/blbx9ak8J0pIVWzl2vqc4fanKRaPxclIeFrP/52LXl0h/L/H3/k15qVJpLf13RZeJ+gviTRL04hYt+oVOL/q+8w9PtAaFu/+aR6sOE8W0CkHD2gLHtBxw0gbaX0XMno1qO9T8CQtApJQL9Th2zEXgjIY0JJ4b5EDkkw/WB8kvq/dVOtU5aAxEBJmJSQzBedJNQCmNtI34FfdQCIx4zZd5KIpP8jPblVu/kLhXglsiQnB7cXi+JNnZ93Kur09Op6uWnagz7a7bh/ZWFB2sCJvPT4cuDhoIaFhdUjjB6xtFMm6dULGXkOjVAHD6mzfmyQ0GqyXi1YH9+rGps0aPhtSQmXnSTke60micW1WlyCAUKKRiOeSaX7DUmoJaCecaP/iyQ1lO3w1c9uggGQ7SMj4rSZnRVCnjgRJ3E3abjW7BnfV0I0myLouNjQ0pIU/B0ZAW69Fdi4+JQ7gU9CTUB9k52kXdKrG3rt+L7NGsoX5XG9kDGJhKEBRe+fAgZKwpBzgaPzyIj7/0kM/7npnElf/fRJo1/aNtPhCJ5DkzUkSTX8fsH3YlHc/4uLIhHn5sRh49/rev9LPkM/G4a5rD7abeD0ablHJhqUSnDxFb+B+j10sk7qaCcChojRydDvROJuamLoOF+F0ucKOZn0/SaNvmvEZSUJeX+NhhCsVotXA9Pufm7z/xuqjpqUvkNFd/5iUSQuZ+Todug4daf+EUK1Ku0cG3MFdn2nSDdp2Ct8Avpt8qUv4+wzM+I8mpsDMF5TxiGSCaiJ5y+ZvVbpp8+b9KB5TV8i9/qnhNqkyc1r+yq1v0+na6wRlxUJAXefXKlKhwWYRrawsHqmDCUYn5t2Pmi1kqAzpVIRtdGPQdPDyONDr07kBKSts7MursgZF/2gl9kVuj2Nhrs3PovRURc60f2I9WV27QJw5kxyB/c7INPTQrU4+yWfjyRJo6WRblun83T6ro1+rWLq8/s2B+2XTilIa8BlQ0L/uetkCN4zA+SVipOW/J12HAdz37uqnS2sPM2loENOMkpNPTj7klVvS7ontttXZ/3BuBO6EdF3CC4uymcSn7mv9Tpw/Ljbv1IB7rsP2Db7dXdw6GHo7ZqA3aTgehAiEV/+6Kcfpk4m0O3oJM1CaU16VOb5fPsmJQychElBZz1hm1Wlt28X6aLrCgFxW29mJp4Az+N1KT9KgHrdXc8fxEPOF/2fJUlB319Aj652Evl93T8+1H+TiBh6drT1WJFibs49w9275bns2CFSelvlWWCxFDeWfTWMDdXk40sjDZU0aXvSqKVHNn+0JPp1puhjtcrDjqEzh1LAwEkY6kQMMuuk5FrNzWTQaqV+PhysdBI899FOF+aE6tCTr/pTGoYGT/2/EEn/79KSk7o8rlsmVadnlSQR2TdYFsZvQ6EgavzkJPCqV8k+d94J/JuPbcLk5CbceiuwdfHvXEkL/8Z8dRNwTgz+vhYC9kKMkN3G7aHRUI92bLcPPRr676F9+XvI6F4n+iKhMeb5AD4D4HmQAk/3Wmt/zxjzIQA/B+DpaNdfsdb++VoaRKLwM0HikDy+dOKzHh9Pzj6ircjJ4OxDfn/pxZzohYCAIz0lL9Xl9WprQHwA0wkKPvSgsrAgr7vvBj73OeALX5Clw5aWgH37rsePjH4v7inVDyUkZZJI142AnaRe0jZNKjpptGroq6gcxZMM8SQiJqETUdeBfiVhC8A/s9YeNsZsAvC4MeYvo98+aq39SD8n84P1JAntPr2+XrnsHCh+Z+Nn30fgl8D0QxXc1omAGkmDe4iQ3Mb74UwM5pOGsrB6QZI05MBVr7t8VY1CwS0IescdwGOPAV/5igwKDzwgrze+EfjIR7bh+sJi3DPlw+/sWtr0S75eyeg/WE1E32DXv+uRuJMBHxrVe72HdaIvElprvwfge9HnZ40x30RCZbU1N6jgZrIzs2N8XGwY38HiP09t1ui0Nv3/accKsLrP+M+XQW5dehNYXd/I98YS3M77mJ8Ph+Q0upHSV+F12hxDOiFzjWl/73mPsxc1/vRPgbe/Hbj+5vrqoGYnFbObVzTphjo9eH9baNRtt+PZMewAVHF0mlKStzXp3EltughYs01ojNkBmWX/VQCvBPAeY8zbAByCSMsfBI4JljzU5OBMBtqBO3eKBGRRJa15AKufk47F+XE5/3r8HPo/2m2RWidPxhcjYufVWTlcw0GXodfXqddFAs3Pu2lP2vTqVTiEkMs59VqbSD7Onwd+9Efd4kohx2ClIs/7W7NbsH3HFmw6/XV3Al+yaMM7RNC1dubQwwiNur5tSGnsS8JebVLd7l7t1JSwJhIaYyoA/gTAe621C8aYjwP4dYid+OsAfgfAP/aP0yUPjZGSh1qTqNdFWrBjU41jiYkkbyQQJlKvKnzIucaAtk71YjUHXoseR66J7i9Lxs+1mvQNBshPnZJzNRqrbcR+/lvGUOt1VxVgcTHZ8cN+Wq+7z7mcaBlTU3Ifd98t9/ubvwm87W0vQf7Y30hjKxX3MDmC+Nkk/sV6efAh0gLhEcpXN0OECR3XbWDo9PC7SdEU0DcJo/L3fwLgP1trPw8A1tqn1O9/COAL/Z6XIzmTjFlSglJGP/8Quqn6QPdnRwmnB9V2W4jTaDjSMDzGtnIAaTTE++g70HI52Xd8HLj5ZomLc0Ej3ad78Q0AcVV4bEze2SYmOfigCjsy4trE+OrCAnDggKSzMcTzvvfJ+z95/Yh8GBnB083N2Jr7vohy7Xrmu7a/eu2oSZ6kkIriO4Z0vmgSqfztnYjYaUBhGy6Citqvd9RAKql901r779T2bZG9CAA/BeBYP+eltJiYENJx2g1VUEomXxX1/QMaSWpo0u/6f200JOH6+HFJ62L5lbk5p85RWi8sSBv5W7MJ7N272mlHm6zdljXUjxxx12aupy6f6PcXfa8kD+vd8F4oDTv1k8XFeJJ6qyVt/9rXZH3306dd2x59FNi371q88sB5fPZPhzE/D7z5zVuwuVZb7ezoRL5+Om9I5dSGvCY94S95pt/9ayd1lBAJ12Mn9IF+JeErAfw0gK8bY45E234FwFuMMfsh6uhpAO/q9YSUPjoEwYJKesUtYLV5kGRPh65BdNJS2m0Xgz51SryHNH8aDZESVE110V96chsNUe0Y29TXIpG1fUkpWqm4DDDtWNIFmwDndCHm512h7FDMM4SdO91gobXLdlsISAfY0pIMQvfdB5w6NYyvfEWu/cY3QkYTfSOaKEByXK4XJNkX/suf0hRyCiWNYjx/EgEvkUOG6Nc7+giAUDWYNcUEAXkuupwg5xDqDsjOqeunJJkS/vNM0nZ4LG29+XnRso4fFwIePizSUBNDr7ilZ1gwk4sOEjo+tOeU7+WydPTFRbleobC6uh7bTMePD9p+XC6C19IV4JL60vHj8nw5w0M/k8ceE6l4003yXi4Dhw4Bx465gfH0aWDr9oLLHfT/TD/dy/9TekEohzVp1O2F3Ekq61qItxbjvQsGnjHD501puLIgCeJOL75rScB9/O9JWoX/fzUaYp/NzYmUO3lSiDc15QjYzdTQnt2lJTmOmStUoXfvXp0coPM65+bC5/TXuSBZqYLyWfmeeH1cEhoNV32b16b0m5qS4ysVJ/SoCR45AuzadS02V9Uo4IcBdKqRHjW7dXZf504aQQcgrWJYa8pTAgZOQj7nsTHXMbXqH5IQREidDNni3Edvp6dzelokwIkTQkKqeElrXZJIbBfgAuHcfuKE2JLttpOOtNf8WqW8b9/ZyBKKJEejIW2lPUq7rlZzUlELjCRpyLbyGtrby+dNlZue27vvFqeSlvyb6Q4+e9Zl0+use/1n9GtbhcinHxCPD8356wZfzdUdq9M5LoIEJAZOQkpA7QHlIEgicqRP8o52G1z1e7vtyhHOzMjIf/SojP6zs6srn3e7Hm092nanT8clOIP9bDs9qoyFzs/HVUsSrd12eaDttvT1UJlNekgZQQCcxsA26GfB7Xq5Bf+c5bJr+9IS8PnPSxB/clLOcf/9AHANDhy4Bj+6pyIPcmYm7iH1SagfSpINl/SH+oH3ENF7Rcjr2una/nEc5Z5L6igL7vqqGMGOw/83ZCP58J8rk5dpf9H2O3pU7J3ZWenkukN2+p/87ezQek4sr0uniw5bnDolUvfYMem7nMCsHVHf+lb8Whs2hP/3VstNBaTU1YOY7i/8TjL7aZW5nBscRkZkH1aKe/BB8eoCzn6fmgJqtc140U7VcEpFbcz72fSdApm6o+sOr8/XCwE6ucJ7Qej4Xj2BfWLgJKR30a+QRviaQpJ9F1Ib6eU7c0aItrgoquKRIyJxSDzaZDzOV32TBj7dN/yK8CQgZ32USu76X/mK9NX5eSe1SLp8tEy8nyN67lzn52iMG2gIv34v70Or+Bpse6UiUq9SkcGCkvaxx+Q767K+6lVCxKXbh1GvvwgLC8BLJr4nJ240nMqhRbS2EZOg/+wkD2e/0k8f0+36fmypm+ReJwZKQtoYOsuD9pLvngecJpM0CPmaCiDnXlwEvvQl8fQ98YRz6mmnhy/59P/USWtinjMdJtp5VCi4at5zc+JxPXJE9qlWge9/P34+Pye0n1WGOS9TCxk+r5Azq9OgwkVv6nUpAkVP7syMm4VRLsszfewxqeJdKMi+07u3YdcuYGwCGK7XRT/X5OtEgJB9xm2aGPqPWyuS2uFP0VmLTdsnBkrCdtvN/9Nl5kMzJZLW39Mk4oujPCXesWOiTjFxupdBrZeBmp+1t1JPtaL9deqUtOXoURkAqlXg7/7OnaPXBXDyeWB5Ofn3pGXT9MpRjLN34sLUlNt/xw43INIcaLUkn5eTojlbnyU0JiYYv9yCrXvKIj6p7hw+HI+NaIRiT/qdFycRNSF7gd9xQqO5Hw7hNS4iBq6Ozs66ArTVqiORDtBTfdKdQQ+YurgunQn33SezAubmXCytk0qb5A3vBbw2JToXEOW0omPHpKOSDAsLTu3sFcZIG4eG5Py852efje/nS0x/9SadBkoTSyelc5+zZ+W/AaQGjZaQ5TKwb5+zdycmhJTNpnhRjx939zY2tgHj4y8BADQXgRtuuknEJ7MM2DCi259Ct68mRzep2M2+63S9TttTwsBJyGfC4LC2S/QsFW3P8LtO5tf/49Gjrpq770jzr91Jte33PpgCRhuqWHTxR+0AYYctFFyYLWkxJELXUeVgtLQEvPCFq504GiQlydhquYRzTnliKmi5LNoCnylLjDB/d3zctZ2DH5PTi0UJZTz4IPD7vy/XGh2V/YtFUVVvvBH4zuIWXLdnj4xM/ggaegAhm1D/mT4Zk8DjfBe7r9boz0me1JRJOTASvuxlwOOPy/ObmnLSg4kYdNNTZfXNCUo8v+oa1SJWZOu0ehaQnnpPErKUBGN/xaJ0vieecNOZdu+W9LFiUUymU6fErsrlRN30VUref7XqzKM9e+T7o4+KdGQ/9EnH71pCPv/5oqaXSsAttzgu0HM6PS3b9u8X6TcyIonnuZzLkfUT2YtF4OGHgY99zAm548flvVIR5xi1hPpt12PT6IwLhIZGSF8l1Nt0R+Co3Cm7P0nF5W9JNl+S6vRcISHgbJhm0xWk5ex5mg1+0jbgOg0DzdosaDScq3929qI5tBKh45mlkkjFsTHgda9zi47Sa9toSKekVGw2gVe/2qmFPAdXgKJJtLQk5NizR7yY998vz6/ZDDt3fBX1u9+V9xtuAB56SJ75wYNyztFRJ53ZVvpW6GSi+QBImx54QEh28KAMLsePr64o0WwKSffvl9cm6rAMjmoRrBGy0fzf/ETdXtWbXh0DocEgRQyUhFpFmpuT/4K5lVRN/dQtwM3zo3ahbUhmfMzMrLbz10LIXsrQ8z6034DXnJwUm2rXLum8TOKenxdybd8unfL220WNrtUkME6V88YbZR/eH8MquZx8npkRycrkg6TK2yH87d8CL3mJHPfFL8p1+LzHx6Wdo6POVmepyLk50V4OH5b/bO9eOW7nTuAXfgH4xCfEHmcSQ60m9zs1JeeamgIWattQKMhzyfuSjC5tjmg+EYHVbnPCl5pJ0CM3pWEvTp6LQMSB24QaCwtu9jm9pJQYOgbnF37mZw6IzDQhQYluIaa1rpDr9xE6kW6+WTrmCyeXI924hQ0lAOUcsKOC802D4dwFYGEBr7yxjOW3bQAA5BvPOjHfbgNLhZVE0a27c1guDKPZlHOzxH6pJJJmerq/tn9dTZ5/yUtcrP3RR2UAGB2VQXHHDpf0MDUl4R5mGOVy0hZOR7vrLjnu/vtdsJ8zU44fl+JSrJiwcyewd+8NOHjwBoyOAhvmnoyCj0vugWoi8pnoh97JoZKkgvrv/uckB8JFcNJcFuoo4IjlO2EA57HTtYfa7bjdx32Z26jzT7l/p3asFZs3x8tu7N0r19qxA/gHd9poIqKaJsEbmJ/HsD5Ro4G8/k7pUC5jubxJBEPD9UcOWIWCTC86cEAKON1/v8TvtHDodXEZEtIYeY6HD4v0oxp5663Cj1OnnDYCCCHvuUeqtl03YXHdBFCpGJw9K47QkPOSduXhw6LO3nyzlGDcu/daXL+zELcX+Tz0u09ETSJNFG3P+E4dn3A64babZByUJOxQ8rAO4LMAdkDmE74pVGOmE5jxNDsrqhC3AfHnwAyYM2fiz4tOmlzOef9OnOinBf1j0yb3n3O2BBOkX/c6iFjyKyrpTqVtGt4wb3pkBOfaw5idBRam4n2LEp6OD5a42LNH1MG9e6Vj+/fv24udSMnxghrFI4+Il/f220Wt3rPHlZdkZlCjgRVX64/fPoajRw2mplxmEBPP6VhaWpLJxHNzkp+6MinjxufhhoMHhd0+GXzi9CK1/BHZz0HVCKlLIcnbS/5kj+hXErYQLnn4dgBfttZ+2BjzAQAfAPDLvZxQS0NOnD192iU66+fPvE/aPvSgUk3ivkxqHh2VkftiOWe0YOO69q2WqHHPKz8rEjApHsXeSND+ifz9P2gM49QpF6ujOr646GxKam2sTMdQwk03CUEefljszFOnws8gKfgfIufmzTLw3X+/3F+1KpJx+3a5Ha7y9M3ZrXjR5AWg0cDrX78R4+POPGD6HnO9H33U3U+5LPdDB1Pp1i14/u7dbqTplhTOSZ18pjRkGdjUcS2qVkCYwL509O1MPeM6BfRFwg4lD98A4NZot08DeAg9klDOJURst+OzCjgSU1UFXAfkM+ez0hX6qFX46VppYuNG53SgJIhds9PKL5Ho+mHDrAwk1SqwAeeAYhHLyGPutHRKvQY9B5u5Oee30AM+VVT2s337XN1Wejg1OmXf+HjmGZH8S0tyrj17XCL6yZPymRLuppuGsLQ0hLNnpZ0TE2L77dwpgwNDGMWiEFGbD2fPOk/56Ouvw/Dp0053JXwS0lWuRygaqLWai1PpOWXsUFrykdAcVULQZfZSwpptQq/k4fNUjZlpiLraF3Tu48mTroaLziVlTJF1QPncZmedfcJOubDgZq73k9nUKyjM6OFsNuXzHXfIq2PV5+iP53/OgbyBDViYXp3TqgdwwD2bWk067dSUy+ucmXGXqtfdQOGX9Ke0e8UrXNEp2phTU8CFC+48L36x9G+GW6pVUb3PnhVpy9We2m3ntAGkLSQ+JyCPjUls8sABIeXtt8v/tH+/SxZot+Uc9TrwmtG6qPWc5Mk/Wy/XzPAGcyCZ+cGKWpReetY4Ryc2kJqJLtHO/8v//1Iuhb8mEgZKHq78Zq21xpigtZFUd9QdKypSoyF/JmNjuqAtBzFdXYzShOohCUjPXNrgopwjI64cY6Egnevuu4Et7aeBpVbyCaKiNeXyxhWpz7bSEwzI/avdY76d+Xnpm6dOufmKeslvQPZbXHSL6XBQardFku/d68hAdfKJJ8ROe+EL3Tmp3rI/U0h88YvyvV4HnnxStuXz8dxTgtKO8xP/+I/le60mDpm77hJn1uc/7zhVKkUnOXlSmH7mjNzQ/LzrBLqQaqXiXrwZv2QDCckZ0UDcYNWLSfqFfnzbMCX0TcJQyUMAT7HimjFmG4CZ0LG67uiBAweCRCVpuM762JgriESiUfugSqqnrbHCdWj9whAh1+MZ5TxAFs29+WZgS90C063VOXfanojYMj/n4qP+/8v+wOyfkyfddCmSlXVttHeW04z4Gh+X/W+6CXjTm0RD21w8JwcfOeIWUZyZkaDme/fhmcbQivZHAcHaO/Pz0rf1LKVcTmzGXK53FXfTJnlECwuyJgYF2S23iBpdqwE3zH0V+P3PiH4KuD+X0zg42rChev0EPmfWzmStR3YgTgDlCMdJn5Sy4+Pxept+0HmA3tFgyUMADwD4GQAfjt7/bK0N4r0uLbmRvtl0Kh9d45QgVEtZE5S2ku9ASxNsIwv/csrPKw7auNvWz+4olXCuvEX60bSb/8f2+QP1E09I5z9xQp7DzIwMSgcPCl9273Z1R7lWx/XjF1a5ii8UNmBo7inRHY9FkfjxcfGwMMgH4LtzGzFz1C1g02xKG6an5dqzsyLlzpyR01M9BcQb28kM9qETz/N54FOfknv+4hdFQt9zD3DDjaPCyLvucqvZPPKIW0/9zBlXFIcLbezdKyctFFyKEh+sLk3H83FWByecUpWl6uWvq0e1bIA2YVLJww8D+Jwx5h0AvgPgTWttEHMnaV9MT0uH2LvXFSbigir0kOr1CLWXvxv51ioF6XijCpfLRWEVig+dNkNE6s7ZsyKAOC+PUo81SWm60PvJ6Vcc1O+4Q7izbx+wtXJOOk+xKCd9eNqVKeBoVKlgiA+DuiSXiCoWcT63YcXcomOHledaLSE/twNuniHtQL5PToq0Xgv4ny8uyjnOnJGBZ+H912N0/F1oRtrAiRPAHa/9Sdx8MzC09KyMEMy745wrVcb96Vmz4l0uFIDihAupsKL7hnpdAqv0enEA5YjOP5ojo142OiX06x19BOGShwDwmvU3h9eJB/FnZyXoS+2Akk4H6v3gfi/XWCs444DX271bOiFmF50o81WXQgFPLWzA8ePidDh+XLSsYlHsst273dLdPPTmm11nqdWAa0rn5KByGWgVgWNnpONQdHH1VHoCyWpdtoAPsNXCcmUzTp90KYDz866WDcEBj/e6fbsMAqdOiWfz9GlXlAqQZHLt1OkV+j9vt4WMv/ALbmHY0VGXnvfgg8CJE5swOflyTL725ThzBpg55Z4b0xYXFlwhY/5nu3aJ3b5nD7BtdNnFdNiJDh1yz46ZH/R2MXGW0iAlrNk7erGh/xTOjKC7nv1ck64flXOtBGQyNHNd5+fdwDu09KxrDBsYNfZ87Xk4c0Y6LkvNs7Q+E83n5mQw379fTlOrARvmv+cewNmF+M1T/6zXXQ/L5ZynihnX4+N4dmlopaocvchcpIbPkEkQDz3kuEuhSvX/2DEXAsnl4iX1tQe32/LeSfCPGR52avjUlHDi6FFpJ7mga9HSrGP8lPFkHfKbmRECPvAAMDubRy73CoyOynO/5RbgOq5fd+qUm4HNupRnz7oqXf3o3l1w2ZIQiP8pepQEkmOtnY5dD3iOfH61kJmbA84VNmFDQXkqqlU8tbgRCwvA2SfcuhNse6HgJsNOTMgof8stwDWNJ6PCO6N4sr1tJTSwufR0vLfNzmJ5xwswMyNtuKb1NM5VtgJ12eX0aWB6Crj/34kpyJzQ+XnX/kpF1Nr9+8UsfOghVzQql3OhIKqCjYZIOU6d0nFHVhHo9Vny/xmOcveoyRQKcUnKWPHCggwiTND3qwP04gOgYPsv/8UlwnMZg717hXe33/5qvHLi684Ly2fOeBSTm1PEZU1CDT1Prl8v53pUTx865jsyIk6Sm24CNhQuAKUSLtS2rtirtJe05jI6Kn/4Lbe4+YFjY0KGzcVzQG07nl00K95fICIvd6bTpdlc6ZytFpCrbsWZU25GfKPh1DGG2fxiUVu3Ssc+fFjaRZORg4VuNwtVDQ9Lf3zmGTcwFQrAD6Ikxeuvl3t+8YuBb3wjfj0973FoyElQPTUq9Jv2mLNN58+H/5/QgMswCgtTcRDU4a5HH5UBCABe+d4J+VMWFuSBMualR56rRRKGoNXUNMnVz3UBcZDs2SPvzHU9X9yE2Zl45kq9LvYizbKFBSHtdaWn5I/lyp3zrZUKwrncphV1sFgENpWXgVwRy22DfLG9suJoMXKInD0rMxNYO5UmTKEgKu7Ro9KRh4fjnffpp93nF7zAhdQ44Zhhj3ZbBhsS/IknXLkNascbN8q+xaKr2r1xo/z2zDPy2+SkEPj0aRduqdXihYVZBZySenJS1GDAqbxJeOlL42EeLsw6MiKaBqUnpZ+usAHI79PTwP/1a5tRrQJ7927Fbbe9AM8787h7OLRDkurkrAFXHAmBS0s+Qo/6XPSFTpNiEbiAoZhPhOrZ6CiwtfxD6b1npvG8Wg1oFYDpRenVp07FJ+q129hYrWLjSrpQGU9Ob8DMDNXI4RXVko6pqSk5lQ7pMO7GeCNzyX1JodMqR0aEwLS9afORLDt3igRnrVSqrTMzcp8M3d14o1zrnnucVH35y+X4e+6Rtj70kNPq6O8YG3NlNXi/k5NybmbzPPywHDM05O6B7efSc3v2yPvx4yLh5udlkJqYcMkQjYbLe2Wyw+Ki7EdyMkf2J27ZJeoCyz/46s06cUWScJDI5aTT79rlaq4wG4SJF+y0W2rLMuyfjjwb7LWRC/387h9Buw1saD0rPbNaxfLoNllI9KQLxh875mJ1ujIdS+DrOXuA8J3mDB2mDLr7SQxMdF9YiM/op73FJjO0Vq+LN3diwlUFoFeV6nOpJJLzne8Efvd3nZR961vdLH6W/ODka0DOl8vJbyQ5BwmGPffvd0KIPqjpadmfqv7oqJz3/vulbZxVw0wsxud5n7WaS/rQGXGPPSZteBabsGly0iUKMMyTEjIS9gmGkLiYDOBiunxtyJ2XjVNzroZFFNl+cmYIZ486QgFAtboJpdILceJRV+mBjgp2CnYc7XxlkSWd7M6sIt1RSSCqt+yArGIwOSn7MyuM07FaLZejSmHN+2WnZUYY4K7JMhhc6q7VAn72Z0USjozI7y+YuIC/9yrg2aWhlUGB4SZ/jRmdgLFvX9wu530D0m7mezebEuOfnXVhxAcfdMvbMf2PAwBj8xpLS+LEAYBbb70Ww5w9vX37al12HchI2CcoXTibgWUqymVgQ3FZehDrLLZaIjIqFXx/Po+jjwjJ+NJJ2bQhgfhq1MoPEytalcuJ9CHxSQZqsWwr33X5lnZbxoQbb3ThxVbLZXgVCk6aAm4yAsMCPL+u/sbl7UolII9l/HApv6IWckA5eHDFcQwsihjbRFHebDrpUioBxeiBFwrAWNmJREbfKaqpAkQG30tvLLqMgmIRP1zKY27OVWFnkgRPoW1RersBJ4kffVQucfgw8MF3bnfX54iUAjIS9gna59r2Yf9YcVcyjWdsDOdLmzE7Lc6Rhx+OT9MCXByOzjctEXTFgJAJQm+mDg+wP+rZI34OLeBmWACuqFO97qQOpa5OFuDCrbz/XA7I5yRwujEizfnWEGbn83jwQWnLHXeI9GOSeLkM5FvnXe+n15dGbqsVz1xhoylSKcbZMD0PUz+oVguo1bCxVEK7NrxiC9brYkpotZwqLf8PSlrA/Z3tNlbP0k8JGQnXgHZbvG20t1aymOYb+EH5WpydvxbHTwKnv+QCx0C8Ng4QD3LTcUIS8p3wS3WQnJzBoR01VCd5bo78XBJtZESkG0d9PySRy7lAPh07DIzrmRqFAla5K+fmxNbT/ZQLyADAcNECrZwzNvWIoeen8Sb5rnVUnbdJ8eXP1YpUg+XCMJqLzlM9MRGfBM7/pVx21fm0kNMraK2sd5eRcLDI5VyBIhLj7FlZ5OXzn9+CJ56IZ4oVi66KmZY8gCt+RFVU9zmtquoZNaH/XsdNy2XpaLR5tMOGdtT27Y68PC/tJG1zUs1mriXJR4noP5hnFvMrGS4UaroT53KIX4A3xDS/cjlOQP3SycH6OH7mg+H26II6l56DFSW85jsnXzB3lnF52u533AF8ee4VmJgAbsh9283XSgEZCXsA3fp02x84IH8O06eYQ8k0Tibl0xHKCtzaCdBuu/xG3dcIdhx+9nPCQ591jR1dzRxwml9oPqueTM52Uu1m1paekscYYr59Acu5oRVJy/Q7PSjQOUOP7srIoBumb1q//OpeVFs16fzyE1pStdvIFdw9a0eWr8ZzsKJ9e/So/DcjI+LQefRROeWddwK//dsvwEbahikgI2EX6LjayIgLfrPcH137fLFCNeA8qTrHVfefbql3ervue2qgX8lkoReUthqFCrfp5H89l1Wnf3F/vtNpw7UTtUeyUADQdtUkmM1F840ZOGyPnoy+HNWV8+87l8sjX7CrjdmkfDQ9Svk6fjRy5XM5lMtDaDTkfrZvdzFHHqbz3klGhn/ItWbT5VUUi3AZGikgI2EHaAJyBOU0JGD1wBvK5qAfYT3QfVBPhuCL5pLu7L52ptVZSgEdBtDb/X2pqumpdPmcxTOLUkeGCSQM8pO83Jc+FABYbstD9SU/USwaDOnsfN641iu1JCW4v7YNowfXbjsHbLXqnFJas223nd3K2KIe7FhFodmUUv933nnt2v7MADIS9gHG0HRhKiA+EAOdM3rSSCZnahg1Nh3s1gQg+bQ6qSedr3h1EScZJQOdPLXaisdf0ugqFueWzEqIRpcbYafVJKQ6TfIut02i11YGE+NqsIakoXYh82CONloSRp/bLVe0mFk61EI4QGjnai4nNr8/jYthmhtvFI0oLfRFQmPMfwTwOgAz1tq90bYPAfg5AMxE/BVr7Z+n18TB4WIkhfsLtnQC+5IPEolLk1F92rXLOVKobrED6f6piadn82uCMiJAEtLbOpRbBpCLFUAjmXUck8cCq4kWghZuEvZohyWgXzZBPxB9keimlpGPhX44WOjTcJBirJDJ+a2WxAepxk5MuOUMNrS8NenWgX4l4acAfAxSAFjjo9baj6TSoqsEoQU9ffhmDuBUJKatUfpwWpQvFHS5SE1KHffT/ZYvOnfKZWC4sBxjyTMLZmVGA9ujBwTtKNLkzmMZaLaQLxSQy+VXkXLlPn3vqP9Q9Haqrto+VCfTmUacl7oS94OL846MuBBOLif7HDvmNIjZWZnOtXOnSP8LpU2d/7w+0BcJrbUPR6UOM6SAXqUiiUZ7TwfMSaZKRfYN1Zfidh0yGBlxBbQ0GTTJSyXG9ZzUOdfMr0jBdnu1l5Ft1HHOFaKr8MNwMQcWadChmFwOThf244bcWcdvtGdJi9/oxshXDlosLcNT6qJg4+PumTSbbsI9Bxx6xHfscLWn0kBaNuF7jDFvA3AIUqE7WAK/W8nDqxXdpCKJQu8mCUhJpj197J/aGaPLT7A+Dm00vxC4dkYUi4hLolwuRlqek21iW4G4cOIxFzAEFIYw1D4PNJsYLnjdr90Glprx7Bn98m9Se8Qo/jVBozghb4HOFRaqInTZTN57uy1kW1yUanAs5fPoo05dTwsdtPSe8XEALwCwH1Kd+3eSdrTW3mutPWCtPbB169YULn31gE4RVnfTaWfaicL9dCaXBj2pvtOGzpdrahabKhYbyxb59upiMe123ASjCqq9wFp4ab6wnRdyw67hOhCvS+bxxLqQkM4c98MWvoT0jGnNVTqvdNVEDlS8D9qFhYJ753lyuVTzt9cvCa21T/GzMeYPAXxhvefM4MABnoH02VnXeThiM8BMlQqIZ9/kck5ScX9OYcrlgKGCdT+2Ee/AWldtt5HL5WMSmaqmrvdCLVInA2i025Cb0swNlczzJSAlpG//acazzcUiLrRMzOHDgYleXUp73kPIcTUxIc+U90S7+7IiIYv+Rl9/CsCx9Z4zgwNtGUoS2igsl8j8TtpmzGXVyxsCLuBeKrmyo0MFGyYekOjKpLrGynd+aAKIz9jgqYYKFs2mixEWCnnkSZymp4ICjlQ+MbXjxZeCWmcvFtH0hKaOofqhRj/OSrD4FqetsUn8ngb6DVH8V8jCLyPGmDMA/iWAW40x+yFLpZ0G8K70mnd1oBd7kFX26ImkLUcXu04N44iv7UMuhT0xIS8Acfd/p/iBQh7LKJXyK5k4/vQqwKm8RKUiGzcWi1jOyfzBfOt8fP0CqpqA0/10TMEPS/ifdbAvaozmKAcsnQ1DaEcrT6HzTB9+2KUXjo7Gs6LSQL/e0bcENn8ypbZkSAD7H3MfuVwCA/W1mozYtOt0BW12REpAzg1cQY/k0xjKLWOoKA3bWARQKOB806xwuVKJq6TlMoCc5NflczmRgBxZcrn40spa5GvJCKwOPGpoTxSwoopqAUpzk5fT9urOnc7O1jFVFkJmCQwWe76s1NEM60MvUpBhCaqAdCZoLY2SSU++pQeUatbYWLRWRpLNF4GpZUAUOCe8zk+p1lLmHOczc81BqsatVh71+uaVoL928JTrmzBMA7XddpntZIRWN3VtCv2Q2L7IyNPrlGghSrtVl7egtqDr7fB0eoGe0VE3IyaNdEQiI+FlDHYMPRNBp5axj/pTjrR3lE4cOlGW2wb5qJctt43YgxFCAoaxvHzOxoypCy2zYp+y9AonJ9PHoqG1RWbhANqhk4+8tXnUattQGdmGLbtsPGbIjPlcLr4ajfa8RN99bZWxQoYodPiEz0w7wXh6Xi4qkLDyvEN5r2tFRsLLFDooTxWTnYC2n65wRrtQB/VJ2FVxxJyTdGHihcBj8mg1XR2bmRmZqsUwBWdUMDsFcETjIj/af0LVWk+9Yl2bWs2gXB5GoTCMalUk+VC97taFY9k5n3HtNgqF/KqBoN127dS81bmxep7k9LQMLKx0KG1yKXxpISPhZQp2iFLJLczZajlVk78xbsh1Ev0AuSYlEfJv9OIc1ebazIyTfFTNdJSBFa4ZUajXXZ0XPX+Ri6xSjdYSnIMMwwKyONAwNnJkouHL1UjpPSmVUKpcE5uqRXCgoEPLrxygwyyzs7JKFGv5AC7G2Pvg1R0ZCS9TUH1qtyX7nx0GiKugJKBOQaO2RkGhpUyIjP58RyDuyqcayf3OnJGOyb7fasm1uGAp15bU9UmZs6knDVOaNxriWKIqqJ0jlFR0SrXbwHJ5I/KjUaPoZaH4jaLuQ6USqtUNMUcrs1z01C+9QhPvfWlJ0tIefVSkvH6uSUkQ60FGwssUeiGl+Xnp+Mzi2LXLBZHprCG0U4HkaLedo8Z3MOqgPrf7gWv6OwAhmXbRczmz06ddZW1mn2jHiA7rAc6/QhJzXiTvkbWQWXtn3z53f5UK8PzqkjspxZpmyMwM8tUqNlWrQDSTggWFOUDt3CnkZ5ErClNWGT96NO7c8lX9tJCR8DKDHv397Yx1sfRlrRbORmGn13ME9fFEKDHFbwfPw/PSAUO7an5eltdmzVJdcMpvlw+ej+FBPfGWa1PmckKUQ4fkGiMjUie1UNgMANizZxP27NmKbeVnZEFUDK1UPSg0gNaCm0e4tCSfNfFHR129Ga4fcvasmzvKfbZvd+X6mV2XFjISDhCh8EShIH/26KgjTKkkHW9kRDrDnj2u1IJe6Vmnk+mZ93pmAxCPfev0TZ8oJOHSkisOzEU8mTQwO+tSwfy0TqD7TBG2Y2bGmXhcNJdSa2pKfj91Sq41OSnTigAh6O7dQK22Gc2mPJsbb3TteewxOY6V1FhFrt0WyT066hxac3NynbNn5R5ZuoSDHtVqSvu0kJHwMkKhIH/25KR0cNaN4ZLrWq3zp9DpeBg9jL5t5e9HEurMME1UHjM3J51+agqxeYTlsqhtx46FpV+IgLrN3IdVCgoFGVwoqUZG5BnQS0wpyyK+JD0XkRkbEwIx5PjII64IFyDHHz/uPKELC7ImBu8FcFKOtWW4cCuwei5iWshIOED4U5iqVRnZJyacRGMJesbjdDY/80EpxSgZdaciQjOBfBJS4ulMk+lpkQqcS6fzQrnEO9VIfV/AaukX8rjqhWC5GCkdNnNzohWwrP6tt8bXWaSkoibAsMgXvgDcdpubJUK1nZ5ZVj+njaxfelYF26yr4unYbVrISDhA6FKKhYI4CiYmZERnkdqVUoFwnYjufF34l0WWQqmgPgH17AYdVmg0XPl7kn56WiQJr8HpPprARLeSH90kCAcSDgYMmnPhnfFxt8/oqNt3bk6eHSCEmZ0VFXTXLnefDCv66jKfv04/1XMKtVSkfaxn4aeBjIQDBr19rF+yY4eM/CQgs1/oIOD+VDd1wjYn7Oq4oCaero0EuG2UoI2GCz2cPu3UT9YS5f508IRsv36RdFw+78IedAKRBLOz8RqnfgkPvo4ccSqrrypraBWcWgVVcSBeNR1w9WXTQkbCywDlsluxiAFxXTWbHVBPRwLiYQQ/IM/PJIsvtUJg5z50SOw8HU8jCQuF+HLW/aJXsi4vu5V9GTLQaXJUS4HVkyw42HTK9+4EfS7/3NxGx1AayEg4QFgrHY3SRk+SJYmoZuqRGYgTkN+1t5PbCL8jEVoa8vvJk/Hl2tdKuvWuqnzhghxPe5GhB10XBognZwPOTiZ02l6vCKmtGt0GtH7QFwkTSh7WAXwWwA7IfMI3JdWYybAajDudPOkyRrSXkx2KKqc/QmvVkNKBRPRjfTyO0Oeik2VqKu4BXF6+dM8iCRxUjh9f7dzR6jW1hXrd1RflnEfGNnuBzhzykWbOKNGvJPwUVpc8/ACAL1trP2yM+UD0/ZfTad5zH+fPi/0zPe3W+NOOAl0LVGe56JFY23V6bqxX7WElV1Ors0y9PHFCMkR0hfFBLEueBN5L0qCQz7uKaHNz7rnRmeVD31sohpl2GKIT+iJhQsnDN0Bm2wPApwE8hIyEfYGEO3FCOszu3bKNUlB7NPWoTgcCbSV6Opk2Rmmpy9HX61iZkZDLybm+9jWxAU+dEsdMWh1Qd25j1kbqfmxIxhuT7N+kc/VSAxaIZ8ldbrMonqdqzEwDeF7SjlnJwzCsBTZskM+c7sNZE3TVkxjT00IUZoRom5DeQmbMUM1kihUn2O7Z4zoppSCdMheDgHrbxZSu/rn7sUf7KcYMDNAm7AZrrTXGJN6ytfZeAPcCwIEDBy4jZWfwOHcOuOYaIZhOuqYrvNGQYDm9llwngXbf4cPOlV8sOonHchi1mgtBLC1Jalel4tZHZCCeSIssfuL2pUS/90Ai9ioZ00IaJHyKFdeMMdsAzKRwzqsSP/gBsGmTSKvjx4UUe/cKqU6ckNgdA8d0wmivIZHLxQPM5bKon5SAJ0645IBaTYjox8LSklqDIuBawXvuZ82Q9SINEj4A4GcAfDh6/7MUznnV4tlonRFjRLqdOCETSxm6GB110230AqOA6zj5fNzzyXXX5+ZEAnIV2hMnRAKG3PHrJeClliYXC90cOGmg3xBFqOThhwF8zhjzDgDfAfCmtBt5NYNJw4Cb3sRZBiROLhf3GoY8iBs2CPGOHJFjmfLG5GufyGnAJ+Ll5G3tF34SepoSvl/vaKjkIQC8JoW2ZFDQHZjeuMVFt5hJsSh2ZK/gvkNDQkQ6dHRw/2KQ5EomXhKyWRRXEXQHHhqS92IR+OEP13deXe4vdK0McVxstToj4RWC9eRrhs6jVxrOCBjGpbJpMxJepXiuOE4uJroNTmk9v4uQCZchQ4Z+YOyAdBFjzNMQb+rlhBEAs4NuxBpwpbYbuLLbfp21dt0LbQ6MhJcjjDGHrLUHBt2OfnGlthu4stueFjJ1NEOGASMjYYYMA0ZGwjjuHXQD1ogrtd3Ald32VJDZhBkyDBiZJMyQYcDISJghw4Bx1ZPQGPMjxpj/1xjzdWPM/zDGVNVvHzTGnDTGfMsY8/cH2c4QjDH7jTGPGWOOGGMOGWMORtuNMebfR20/aoy5cdBt9WGM+WzU7iPGmNPGmCPqt8v6uacOa+1V/QLw/wH4sejzPwbw69HnPQD+BsAwgOsBfBtAftDt9dr+RQD/IPr8WgAPqc9/AVle9yYAXx10W7vcx+8A+BdXynNP+3XVS0IAuwA8HH3+SwD/MPr8BgD3WWvPW2v/DsBJAAcH0L5OsAAouTcDOBt9fgOAz1jBYwBqUdWDyw7GGAOZg/pfo01XwnNPFRkJgW9A/ngA+EcAnh99vhbAd9V+Z6JtlxPeC+DfGmO+C+AjAD4Ybb8S2k68GsBT1tq/jb5fSW1PBVfFLApjzJcAjAV++lWICvrvjTH/N6RUR4p1tNaPLm1/DYB/aq39E2PMmwB8EsDtl7J9ndCp7dZalkF5C5wUvCpxVZDQWtutY94BAMaYXQB+Itr2JJxUBIDxaNslRae2G2M+A+AXo6//DcAnos+XfdsBwBhTAPBGAC9Tmy+Ltl9KXPXqqDFmNHrPAfg1AH8Q/fQAgDcbY4aNMdcDuAHA1wbTykScBfBj0efbAFClewDA2yIv6U0AnrGuNuzlhNsBPGGtPaO2XQnPPVVcFZKwC95ijPn56PPnAfwRAFhrv2GM+RyA4wBaAH7eWnsZrMwQw88B+L1IoiwhKqwM4M8hHtKTABoA/s/BNK8r3gxPFb1CnnuqyNLWMmQYMK56dTRDhkEjI2GGDANGRsIMGQaMjIQZMgwYGQkzZBgwMhJmyDBgZCTMkGHAyEiYIcOAkZEwQ4YBIyNhhgwDRkbCDBkGjIyEGTIMGBkJM2QYMDISZsgwYGQkzJBhwMhImCHDgJGRMEOGASMjYYYMA0ZGwgwZBoyBFXq60xg7a4x80e+hbf3+xs/6exrn7fe3S3XNhP2tlVfoc6dtafymSxf10o6L8Vvoc1rtwYu+A3xz9v+x1t6JdWJgJJwFcGhoCMjlgEIh/h7alu3T2z7R5+W2QasFtNtY9R7alu3T3z7mP70L9sC9I2lwIVNHM2QYMDISZsgwYAys7qgx5hikYG2GDFcqStbaves9ySArcC9Zaw8M8PoZMqwLxphDaZwnU0czZBgwMhJmyDBgDJKE9w7w2hkypIFU+nC2IEyGDANGpo5myDBgZCTMkGHA6IuExphfMsYciV7HjDHLxph69NsvRtu+YYx5rzrmQ8aYJ9Vxr422/7gx5nFjzNej99sSrvlZdexpY8yRaPtBtf1vjDE/tdaHkCGDD2PMncaYbxljThpjPhD4/X3GmOPGmKPGmC8bY65Tvz1ojJk3xnyhp4tZa9f0AvCTAL4Sfd4L4BiAMiT2+CUAk9FvHwLw/sDxLwWwXR3/ZA/X/B0A/yL6XAZQiD5vAzDD79kre63nBSAP4NsAdgIoAvgbAHu8ff4egHL0+Z8A+Kz67TURP77Qy/XWo46+BW6p4xcB+Kq1tmGtbQH4nwDe2Olga+1fW2vPRl+/AWCDMWY4aX9jjAHwJl5TXQsASgAyD1OGtHAQwElr7SlrbRPAfQDeoHew1v6VtbYRfX0MwLj67csAnu31YmsioTGmDOBOAH8SbToG4NXGmC3Rb68F8Hx1yHsisf0fjTHXBE75DwEcttae73DZVwN4ylr7t6odrzDGfAPA1wG8W5EyQ4b14FoA31Xfz0TbkvAOAH+x1outVRL+JID/ba2dAwBr7TcB/BaALwJ4EMARAMvRvh8H8AIA+wF8D6JSrsAY8+Lo2Hd1uaaWvIiu+1Vr7YsBvBzAB40xpTXeT4YMa4Ix5h4ABwD827WeoysJjTE/rxwg26PNb8ZqQnzSWvsya+0tAH4A4ES0/Slr7bK1tg3gDyGinuceB/DfAbzNWvvtDm0oQNTbz4Z+jwaBRYhtmSHDevEk4prceLQtBmPM7QB+FcDru2hxnbEGo3UzgDkAG73to9H7BIAnANSi79vUPv8UwH3R5xrE4H1jD9e8E8D/9LZdD+eYuQ7AWQAjgzbqs9eV/4I4F09FfYyOmRd7+7wU4ry5IeEct6JHx8xaGvh2Esnb/r8AHI8a/Bq1/T9BbLajAB4gKQH8GoAfQlRXvkjkTwA4oM7xKYjNp6/30xCHzhEAhwHcNeg/L3s9d14Qv8aJiGi/Gm37VxCpB0gE4CnVdx9Qx/4vAE8DOAexJ/9+p2tlaWsZMgwYWcZMhgwDRkbCDBkGjIyEGTIMGBkJM2QYMDISZsgwYGQkzJBhwMhImCHDgPH/A7qcPko5/0VnAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "hgrid.make_plot(show=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "recognized-garbage",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
